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NOTE  ON  THE  DETERMINATION  OF  THE  FORM  OF  AIR  SHOCKS  FROM  THE  DF.CAT  CURVE 


ABSTRACT 

Expressions  connecting  the  first  derivative  of  the  shock  parameters 
just  behind  the  front  with  the  decay  curves  are  given  in  the  form 
convenient  for  comparison  with  experimental  results,  for  plane,  cylindrical 
and  shperical  shocks.  Comparison  of  pressure-time  gage  measurements  on 
"spherical"  shocks  with  the  decay  curve  for  shocks  from  spherical 
Peatolite  charges  is  included.  The  discussion  is  preliminary  to  the 
determination  of  the  form  of  the  shock  from  these  charges  now  being 
undertaken  by  integration  of  the  hydrodynamic  equations  along  the 
Kirkwood-Brinkley  decay  curve. 
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The  determination  of  the  pressure-time  relation  at  a  point  and 
the  pressure-distance  relation  at  a  given  time  from  the  decay  of  the 
peak  pressure  of  an  air  shock  has  been  hampered  by  the  small  range  of 
peak  pressures  over  which  reliable  measurements  were  available*  Recent 
work  in  this  Laboratory  on  spherical  Pentolite  charges  has  confirmed 
the  general  form  of  the  decay  curve  obtained  by  compilation  of  data 
from  a  number  of  sources In  the  following  note  the  slopes  of  the 
pressure-time  and  pressure-distance  curves  just  behind  the  front  are 
computed  as  a  preliminary  to  a  more  complete  study. 

Following  Kirkwood  and  Brinkley^  we  calculate  the  pressure-time 
or  pressure-distance  variation  at  the  shock  front  by  combining  the 
Hugoniot  conditions  with  the  hydrodynamic  equations.*  For  comparison 
with  experiment  the  Eulerian  form  is  more  convenient  than  the  Lagran- 
gian  form  used  by  Stanton  and  Cassen3  and  will  therefore  be  used  from 
the  start,  although  the  relations  derived  here  can  be  obtained  from 
theirs.  To  distinguish  the  variation  of  a  quantity  at  the  front  from 
the  usual  total  derivative,  we  write 
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where  U  is  the  shock  velocity,  r  the  space  variable  and  R  its  value  at 
the  front. 

te +  u  If  *  k  * (u " u)  If  ' u  (§n  ‘  p2  If  )  (u) 

tfiere  u  is  material  velocity  in  the  shock  front,  p  and  pQ  the  densities 
in  and  without  the  shock  front. 


The  hydrodynamic  equations  can  then  be  specialized  at  the  front  in 
the  form: 


Continuity  equation  and  energy  equation. 
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where  p  is  the  pressure,  P  is  the  pressure  at  the  front,  and  ▼  is  the 
material  velocity,  a*  is  the  variation  of  pressure  with  density  at 
constant  entropy  in  the  front,  a  •  0  for  a  plane,  1  for  a  cylindrical, 
gnd  2  for  a  spherical  shock  front. 


Dynamical  equation, 

Du  Po  „  dv  1  dp 
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(3) 


or  combining  (2)  and  (3)  and  utilising  p 
pressure  of  the  undisturbed  air. 
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In  an  appendix  the  notation,  equations,  and  Hugoniot  conditions  are 
•umarlsed. 
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Since  at  very  high  peak  pressures  Dp/DP  is  very  small,  it  is  convenient 
to  express  u  in  terms  of  P  and  p  at  the  front*  The  following  ex¬ 
pressions  have  been  found  useful: 
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Experimental  work  on  spherical  charges  shows  that  decay  is  linear 
both  at  very  low  and  at  very  high  peak  pressures.  The  equation 

(P  -  p  )Z  »  (2.26  +  - -T7T--3  7)  atra  ft  (Q) 

1  +  0.23  lb^/:5ft 

where  Z  *  R/(charge  weight)1/3,  has  been  found  to  fit  data  of  this  Labora¬ 
tory  on  spherical  Pentolite  over  the  range  from  the  charge  surface  to  Z  ■ 
100,  within  the  dispersion  of  measurements.  Calculations  of  peak  pressure 
from  equation  (8)  can  be  relied  upon  at  large  distances  from  a  detonation, 
probably  to  better  than  5%,  as  is  shown  by  the  agreement  between  indepen¬ 
dent  workers .U  At  small  distances  the  non-spherical  form  of  the  shock 
produced  by  initiating  a  spherical  charge  with  ary  existing  detonator  and 
other  so-far  unidentified  factors  introduce  an  uncertainty  of  the  order 
of  10£.  If  the  scaled  distance  is  expressed  in  charge  diameters,  rQ, 

instead  of  in  ft  lb"1/3,  x  -  r/rQ  -  7.52  Z 

(P  “  p J  x  -  (17.0  +  - ~ - --»)  atm  charge  radii  (8a) 

°  1  +  lj.lx  10"3xZ 

Figure  1  is  a  plot  of  the  relevant  parameters  in  these  expressions 
taken  from  the  tables  of  Kirkwood,  Brinkley  and  Richardson.^  These 
tables  are  also  used  in  obtaining  the  peak  pressure  close  to  a  charge 
from  the  measured  shock  velocity.  In  neither  case  can  the  possible 
errors  in  the  equation  of  state  used  in  the  Hugoniot  conditions  intro¬ 
duce  an  error  in  the  results  at  high  pressures  comparable  to  the  un¬ 
certainty  of  measurement,  as  the  relations  between  peak  pressure  and 
shock  velocity  and  between  peak  pressure  and  particle  velocity  are  vary 
insensitive  to  variations  in  density  behind  the  front,  while  the  terms 
in  equations  h-7  containing  other  parameters  are  small.  We  discuss 
pressures  above  and  below  100  atmospheres  separately. 

Within  a  few  charge  radii  from  the  surface  we  have  the  rough  approxi¬ 
mation 
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charge  surface  to  1/3  At  100  atmospheres  pressure,  and  can  be  neglected 
in  coaqparison  with  p/pQ  over  this  range,  while  pQ  is  negligible  ooqpered 

to  P.  Baaulta  of  computations  from  the  expression 
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Properties  of  air  at  the  shock 
front,  from  Kirkwood  -  Brinkley  tables. 
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are  shown  in  Figure  2.  The  "lopes  are  given  in  charge  radii  for  better 
visualization  of  the  shock.  P/^  is  the  distance  in  charge  radii  from 


the  front  to  the  intersection  of  the  tangent  to  the  pressure  distance 
curve  at  the  front  and  the  lina  of  P  •  pQ.  It  is  interesting  that  this 

distance  remains  about  l/$  the  distance  between  the  front  and  the  center 
of  the  charge  throughout  the  range  of  these  measurements.  _  1  $  is 


computed  from  1/P  dp/dr,  using  equation  (1). 


At  pressures  less  than  100  atomospheres  the  perfect  gas  approxima¬ 
tion  can  be  used.  For  a  gas  of  specific  heat  ratio  setting  (B-P0)/p0*y 
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The  second  form  is  obviously  convenient  for  the  use  of  equation  (8). 

Results  of  this  computation  for  air  (y  -  l.U)  are  plotted  in  Figure  2. 
Equation  (7)  has  been  used  to  compare  the  observed  slope  of  the  pres¬ 
sure-time  curve  recorded  oscillographically  with  the  slope  to  be  expected 
from  the  decay  curve  (peak  pressure-distance  relation).  Agreement  has  usual¬ 
ly  been  obtained  dan  the  shock  was  very  long  conparod  to  the  gage  size. 

Good  agreement  cannot  be  expected  with  data  on  shocks  close  to  small 
charges  when  large  gages  are  used  (the  BRL  gage  and  fractional  pound 
spheres,  for  instance)  since  the  gage  gives  an  average  pressure  over  a 
considerable  part  of  the  shock  wave.  A  few  points  obtained  in  this  way 
with  1/U  pound  Pantolite  charges  are  included  in  Figure  2. 

It  ie  obvious  that  the  method  used  to  obtain  equation  (h)  can  be 
iterated,  differentiating  the  hydrodynamic  equation* ,  to  obtain  succes¬ 
sive  derivatives  at  the  front.  These  in  turn  can  be  inserted  in  s  Tay¬ 
lor's  aeries  in  the  distance  from  the  front  and  the  form  of  the  pressure 
time  or  pressure  distance  relation  behind  the  front  obtained  to  the 
discontinuity  in  any  derivative.  This  discontinuity  would  be  expected 
to  be  the  boundary  between  the  explosive  gases  and  air,  at  which  even 
the  first  derivative  of  the  pressure  is  discontinuous ,  unless  a  second 
shock  has  formed.  Although  theoretical  discussion  of  the  form  of  a 
spherical  shock  ie  difficult,”  it  is  clear  that  a  second  spherical 
shock  form  much  more  slowly  than  a  second  plane  shock.  As  experi¬ 
mental  studies  cm  small  charges  have  not  been  carried  out  with  initiation 
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resulting  in  a  spherical  shock,  the  appearance  of  multiple  shocks  In 
pres sure- tine  records  cannot  be  regarded  as  evidence  as  to  their  forma- 
tion  from  a  detonation  having  spherical  symmetry.  The  question  of  the 
region  of  continuity  thus  remains  unanswered  until  further  theoretical 
work  has  been  done.  It  is  doubtful  that  the  experimental  results  are 
sufficiently  accurate  so  that  reliance  can  be  placed  on  higher  deriva¬ 
tives  of  equation  (8)  than  the  first,  in  any  case.  Computations  of  the 
pressure  behind  the  shock  front  as  a  function  of  time  by  the  method  of 
characteristics  are  being  made,  using  the  Kirkwood- Brinkley  peak  pres¬ 
sure  vs  distance  curve?.  Similar  confutations  will  be  made  using  the 
experimental  curve,  if  further  experimental  work  shows  significant  dif¬ 
ferences  between  the  curves.  The  estimation  of  positive  impulse  from 
decay  may  prove  more  accurate  than  its  experimental  determination 
close  to  charges. 


Jane  Dewey 
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APPEND II 


Symbols  used? 

r  *  space  variable,  for  a  spherical  shock  the  distance  from  the 
charge  center 

B(t)  •  the  value  of  r  at  the  shock  front 
t  •  time 

p(r,t)  ■  pressure 

P(B)  •  pressure  at  the  front 
p  ■  density 

.2  . 

S  •  entropy 

U(R)  -  shock  velocity  relative  to  gas  in  front  of  shock 
v(r,t)  ■  material  velocity  relative  to  gas  in  front  of  shock 
u(R)  *  value  of  v  at  the  front 

a  is  defined  V  •  -  |j  *  £ 

V  -  charge  weight 

Z  • 

r  -  (r  -  po)/p0 

y  «  ratio  of  specific  heats 

The  subscript  o  refers  to  the  gas  in  front  of  shock,  which  is  assumed 
to  be  uniform. 

No  subscripts  are  used  for  quantities  just  behind  the  front. 

The  one  dimensional  hydrodynamic  equations  arei  in  Bulerian  co¬ 
ordinates 
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Hugoniot  Conditions 

A*  For  any  material 
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Change  in  internal  energy  of  material  passing  through  the  front  is 
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